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Calculating multiloop integrals using dimensional recurrence relation and 
"D-analyticity 

R.N. Lee a * 

a The Budker Institute of Nuclear Physics 

We review the method of the calculation of multiloop integrals recently suggested in Ref.Q]. A simple method 
of derivation of the dimensional recurrence relation suitable for automatization is given. Some new analytic results 
are given. 



1. Introduction 

Recently, in Ref Q] a method of multiloop in- 
tegrals evaluation based on T> recurrence rela- 
tions [2] and £>-analyticity was suggested (DRA 
method) . In this contribution we give a brief re- 
view of this method. We also provide a simple 
method of derivation of the dimensional recur- 
rence relation well-suited for automatization. 

2. DRA method 

The DRA method has been described in detail 
in Ref. [T|. It consists of the following steps: 

1 . Make sure all master integrals in subtopolo- 
gies are known. If it is not so, start from 
calculating them. 

2. Pass to a suitable master integral j( v \ It 
is convenient to choose a master integral 
which is finite in the basic stripe. For this 
purpose, e.g., increase powers of some mas- 
sive propagators. 

3. Construct the dimensional recurrence rela- 
tion for this master integral. The general 
form of this recurrence is 



fP-*) = C(p)jV>) + R(p), 



(1) 



where C{T>) is some rational function 
and -R(Z?) is a non-homogeneous part 
constructed of the master integrals of 
subtopologics in T> dimensions. 



4. Find a general solution of this recurrence 
relation 



jP) =E- 1 (Z?Mz) + J; 



(z>) 



(2) 
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where Jj h is a specific solution of the in- 
homogeneous equation, E _1 (£>) is the so- 
lution of the homogeneous equation, and 
lo(z) = w(exp(i7r2?)) is arbitrary periodic 
function. 

5. Fix the singularities of w(z) by analysing 
the analytical properties of the master inte- 
grals and summing factor £(X>). 

6. If needed, fix the remaining constants from 
the value of the integral at some space-time 
dimension T>. 

Step 5 is the key point of the DRA approach. 
In order to perform this step one needs to de- 
termine the position and order of the poles of 
£(£>) J^ on suitable vertical stripe of width 2 
in the complex plane of T> {basic stripe). This in- 
formation can be extracted from the parametric 
representation of the integral, either manually or 
(semi-)automatically using the FIESTA code [3]. 
In general, the number and the order of poles es- 
sentially depend on the choice of the master in- 
tegral J^, the summing factor £(£•), and the 
basic stripe. The proper choice may essentially 
simplify the last step of the approach. 

Before we proceed to the example of the appli- 
cation of the DRA method, we would like to de- 
rive formulas convenient for the automatic deriva- 
tion of the dimensional recurrence relation. 



3. Dimensional recurrence relation 



such functions as 



The original derivation of the dimensional re- 
currence relation [2] is based on the parametric 
representation. For the integral without numer- 
ator which can be represented by some graph 
the final formula has the form of some sum over 
the graph trees. For the automatic calculation it 
may be desirable to have the possibility to ob- 
tain the dimensional recurrence relation without 
any reference to the graph and/or to the para- 
metric representation. In this Section we obtain 
the corresponding formulas using the Baikov's ap- 
proach which consists of the "changing of inte- 
gration variables" from loop momenta to scalar 
products (or denominators) [3] . We briefly review 
the derivation of the corresponding transforma- 
tion keeping also Independent factors omitted in 
the original derivation of Ref . [4j . 

Assume that we are interested in the calcula- 
tion of the L-loop integral depending on E lin- 
early independent external momenta p±, . . . ,Pe- 
There are N = L(L + l)/2 + LE scalar products 
depending on the loop momenta If. 



k-q,j] i-l,...,L; j — l,...,K, 



(3) 



where <2i,...,l = Ii,...,l, Ql+i k = Pi,...,e, and 

K = L + E. 

The loop integral has the form 



jvn*)=J d %i ll j(n 1 , 

d v l L . . . d v h 



n LV/2 D ni D n 2 p 



n N ) 
i» ( 4 ) 

N 



where the scalar functions D a are linear polyno- 
mials with respect to Sij. The functions D a are 
assumed to be linearly independent and to form 
a complete basis in the sense that any non-zero 
linear combination of them depends on the loop 
momenta, and any sn~ can be expressed in terms 
of D a . 

The integral J^ (n) can be considered as a 
function of N integer variables. It is convenient 
PQ to introduce the operators Ai, Bi which act on 



{Aif){. . . ,m, . . .) = mf(-.-,rn + l,...) 
(B i f)(...,n i ,...) = f(...,n i -1,...) (5) 

Let us first transform the innermost integral 
J dPli/ir' / 2 in Eq. ((4]). The integrand j depends 
on l\ via the scalar products su (i = 1 . . . L + E). 
Writing l\ = Z ^_ || +h±, where Zi|| is the projection 
of l\ on the hypcrplanc spanned by q-2 . . . qx, we 
obtain 

d v h d v - K+1 l l± _ d K - 1 l n 

~^D/2 ~ 7T (T>-K+l)/2 7T (K-l)/2 



lX (V-K-l )/2 
V{qi,...qK) y d 



( V{q u ...q K ) \ ( 

\ fl V(q2,...qK)J 



T[(V-K+l)/2] 

ds 12 . . . ds 1K 



^K-DfrjpK-iy^ 



qi<) 



(6) 



where V{q\, . . . qn) = det{s,j \i t j=i...K} is a Gram 
determinant constructed on the vectors q\, . . . qx 
and fi = ±1 for the Euclidean/pseudoEuclidean 
case, respectively. Note that 

V(qi,...,q K ) = P(D 1 ,...,D N ) 

is a iiT-degrce polynomial of D a . 

Repeating the same transformation for 
I2, ■ ■ ■ , II, we finally obtain 



J(n) 



n (L~N)/2 fl LV/2-N 




r [(V - K + l)/2,...,(D-E) /2] 

L K M(X>-if-l)/2 






[V (qi,...,q K )} 



[V(pi,...,p B )]P J - B - 1 Va 



i(n) 



(7) 



In order to use this formula in explicit calcu- 
lations, we also need to determine the limits of 
integration over the Sij variables. However, for 
algebraic manipulations we only need to keep in 
mind that the integration by part does not gen- 
erate any surface terms. 

The lowering dimensional recurrence relation is 
immediately obtained by replacing V — > T> + 2 in 
Eq. ([7]) and comparing the resulting expression 
with the original one [Tj. We obtain 



j( c+2 > (n) = 



(2^) L [V (Pi 



,Pe) 



(D-E-L + 1) L 
x(P{B 1 ,...,B N )J^{n). (8) 

In order to obtain the relation between master 
integrals, we have to use IBP reduction for the 
right-hand side of Eq. ((5]). The complexity of 
this reduction strongly depends on the integrals 
appearing in the right-hand side. The lowering 
dimensional recurrence relation © contains inte- 
grals with indices shifted by at most K = L + E 
in comparison with the integral in the left-hand 
aside. 

The raising dimensional recurrence relation is 
more " economic" from this point of view. In the 
original Tarasov's derivation the parametric rep- 
resentation of the loop integral was used. For the 
integral given by some graph, the result is ex- 
pressed in terms of the trees of this graph. How- 
ever, for the automatic derivation of the raising 
recurrence relation this formula may be inconve- 
nient. Therefore, it is desirable to be able to ob- 
tain the raising recurrence relation without any 
reference to the graph. In order to obtain the 
raising recurrence relation we use the identity 



det <^ 2 5 » 



d 



\i.j = l...L 



[V(qi 



,Qk) 



(2a) L V( Pl ,..., PE )[V( qi ,...,q K )Y 



(9) 



The proof of this identity is based on the Carl 
Jacobi theorem about determinants and will 
be presented elsewhere. The raising dimen- 
sional recurrence relation is obtained by replac- 
ing T> -4 V - 2 in Eq. 0, substituting 
V ( Pl , . . . ,p E ) [V (<&, . . . , q K )] iV - K - 3)/2 with the 
derivative and integrating by part. We obtain 



j(p- 



n 



(-M/2) 2 



dPl 



..dPh (2 s «d 

det { — U,j=i...L \j (n) 



(M/2) 



n LV/2 
L 



dsi 
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x (det/2*«|^M. J . =1 ... i }j(^(n). (10) 



Comparing Eq. (|10p with Tarasov's formula 
we obtain for the case of integral corresponding 
to some graph: 



{2 s ^- l ^A k \ l . ]=1 ... L 






where the sum goes over all trees of the graph, 
and i\, . . . ,%l enumerate the chords of the tree. 

4. Example 

Let us demonstrate the application of the 
method on the calculation of the following four- 
loop vacuum integral: 



JW = 




d u kd u ld u rd u p 



Tr 2T >k 2 l 2 r 2 (k + l + ry (p-k-lf + 1 



[( p -ky + i}\p 2 + i][(p+ r y + i} 



This integral has been considered in Rcfs. |5l6j . 
In Rcf. [5] this integral has been evaluated nu- 
merically using the Laporta's defference equation 
method. In Ref. [B] this integral has been con- 
sidered using the dimensional recurrence relation. 
However, in that paper in order to fix the peri- 
odic function parametrizing the homogeneous so- 
lution we had to resort to the Laporta's differ- 
ence equation. Here we present the derivation 
entirely based on the DRA method. This deriva- 
tion serves solely as the illustration of the DRA 
method. The final result for arbitrary T> coincides 
with the result of Ref. [6] . 

1. There are four master integrals in the 
subtopologies: 



■n 



CD) 



j(P) 




(11) 



These integrals are expressed in terms of F- 
functions, see, e.g., Ref. [6]. 

2. The integral has no ultraviolet divergence 
for V < 4. At V = 4 the integral has 
a simple pole. The integral has a simple 
pole also at T> = 2| due to the infrared in- 
frared divergence. This divergence comes 
from the region where fc, I, and r arc small. 
Subtracting from the integrand the quan- 
tity k 2 l 2 r 2 (k + l + r) 2 (p 2 + l) 4 , cor- 
responding to a scalclcss integral, we eas- 
ily establish, that J^ is finite when 
2 <KeV < 2|. We choose the basic stripe 
as5 = {D| ReT> e (2,4]}. 

3. The dimensional recurrence for J^ ' reads 




Figure 1. Pole structure of functions ii, *2 (upper 
row), t3,ti (lower row). We plot \U [T> + 2k)\ for 
fc = 0,...,3. 



AV+2) = 48(3X>-ll)(3X>-7) {v] 
(D - 3) 4 (£> - 2) 2 

+ c (V)j(V) , JP) 7 (P) , JX>) T (X>) , jp) T (V) 



c 2 ^2 



c 3 ^3 



c 4 "4 ) 



„(£>) 



where c\ are some rational functions not 
presented here for brevity (see Ref. [5]). 

4. The summing factor obeys the equation 

£(Z>) 48(3£>-ll)(3X>-7) 



£ (2? + 2) 



(£> - 3) 4 (P - 2) 2 



(12) 



Since J^ has simple poles in the basic 
stripe at T> = 2|,4, we choose the sum- 
ming factor to have zeros at these points. 
Namely, we choose 



E(3?) 



s(^ + f)r 2 (p-3)r(|) 
8 p r(3-D)r 



"iV 
> 2 



2 - 



(13) 



The general solution of the dimensional re- 
currence has the form 



5. The left-hand side of Eq. (fT5j) has no singu- 
larities on S, so the right-hand side should 
also be a holomorphic function. The func- 
tions ti (T> + 2k) have poles at V = 2?i_6 
where 

1 2 

P 1= 2-,P2 = 2-,P 3 = 3, 

2)4 = 3^5=3^6 = 4, 

The pole structure of U (V + 2k) is demon- 
strated in Fig. 1. Note that Si (V) may 
have poles only in the points where the in- 
dividual terms fcj (T> + 2k) are singular. 

Therefore, we obtain 



6 r t 



uj(z) = j2J2 c ^ [ cot 



1=1 r=\ 



■ziP-Vi'. 



-const, 



(17) 



E (V) J™ = w («) - J> ( p ) ( 14 ) 

i=l 

Sl {V) = ^ U (V + 2k) (15) 



k=0 



t l (V) = X(V)c\">J l 



(V) T (V) 



(16) 



where n is the order of the pole at T> = T>i , 
and the coefficients c[ should be chosen so 
as to cancel all singularities in the right- 
hand side of Eq. (fi"5j) . Their determina- 
tion is reduced to the solution of some lin- 
ear system which we do not present here for 
brevity. 



As it is shown in Ref. jT|, an L-loop in- 
tegral is bounded in the limit T> — > ±ioo 
by z ±L ' i \\ogz\ l/ : where v is some irrele- 
vant exponent. Using this fact and the ex- 
plicit form of £ (T>), it is easy to establish 
that E (V) J^ and s, (X>) fall down when 
T> — > ±ioo. Therefore, the constant in Eq. 
(fT7|) should be chosen in such a way that 
ui (z) falls down when z — > 0, oo. The first 
term in Eq. (jTTJ) has different limits when 
V — > ±ioo, and we obtain 



const 



■££■ 

1=1 r=l 



6 I'; 



"*)' 



£<c(+o r 

(18) 



1=1 r=l 



Eqs. ([T5 ]) .([T7 p . and jTHJ) entirely determine 
jC ) for arbitrary 2?. However, we may 
want to find the coefficients c[ in Eq. (fT7|) 
cxplicitcly. Using the fast convergence of 
sums in Sj (2?) and keeping in mind the pos- 
sibility to use the pslq algorithm [7], we 
find that values of all these coefficients are 
compatible with zero, at least, up to 10 -500 . 
Therefore, we conclude that 



I ^ 500 n 
oj{z) = 0, 



.-,1)0 



(19) 



where = denotes the equality checked nu- 
merically with 500 digits. 



6. Our consideration allowed us to fix all con- 
stants within the method. Adopting the 
guess (fT9|) . we obtain 



JV 



-E" 1 ^) $>(£>) 



This result coincides with that of Ref. [6] . Us- 
ing the pslq algorithm, we can express the expan- 
sion near T> = 4 in terms of conventional £- values: 



j(4-2e) 350 e 3 r[-l + e] 4 



5Cs 
/292137T 



1 + e 
7(1 ■ U * 6 



[1 + 2e 3 + 3e 4 + O (e 5 )] 
'As 



378 J V 30 
1820C 2 ,6 - 5038C 3 C 5 ) e 3 



212C 7 e 



V 32400 

^13255C 9 73l7r 4 C 5 2006tt 6 C3 1006(1 \ 



I 3 

o( e 5 )] 



6 



189 



Note that the factor in the first line is cho- 
sen so as to provide the uniform transcendentality 
weight in the rest of expansion. In Ref. [5] the 
first term of the above expansion has been found 
analytically and the remaining terms have been 
found numerically with 40-digit precision. 

Let us also present the result for the expansion 
around T> = 3, which can be important for the 
calculations in hot QCD: 



J( 3 - 



r[-l/2 + ef 
1 + e 



7T~ 

96 



2717T 4 

2880 



- 2 log2-^ 



11C3 
16 

417T 2 

48 



O 



In Ref. [B] the first term of the above expansion 
has been found analytically and the second term 
has been obtained numerically 

5. Conclusion 

We have briefly reviewed the method of cal- 
culation of multiloop integrals based on the V- 
rccurrence and 2?-analyticity. The method ap- 
peares to be powerful enough to deal with the 
most complicated cases. We have also derived 
convenient formulas, Eqs. ([8]) and (|T0| . suitable 
for the automatic derivation of the Tarasov's di- 
mensional recurrence. For a specific four-loop 
master integral we have presented in analytic 
form several terms of the expansion around T> = 4 
and V = 3. 
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